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fSj ' The search for stable heavy exotic hadrons is a promising way to observe new physics pro- 

■ cesses at collider experiments. The discovery potential for such particles can be enhanced 

Q ' or suppressed by their interactions with detector material. This paper describes a model 

^ ] for the interactions in matter of stable hadrons containing an exotic quark of charges ±|e 

' or ±|e using Regge phenomenology and the Quark Gluon String Model. The influence 
of such interactions on searches at the LHC is also discussed. 

^: 

1 Introduction 

CD 

1 i | Searches for slow- moving stable massive particles (SMPs) which interact within a detector 
^ ■ offer a promising means of observing new physics processes at colliders. Their strong 
j> ■ production mechanisms will allow hadronic SMPs to be among the exotica for which 
O . the LHC can open a discovery window for comparatively small amounts of integrated 

i luminosity (~ Ifb^^) [1]. The hierarchy problem suggests the manifestation of hitherto 

I unobserved physics processes at TeV collision energies, and it is thus prudent to consider 

O ■ the possibility of heavy exotic stable quarks. Furthermore, such particles are predictions 

^ of phenomenological implementations of a number of candidate theories which extend the 

O ■ Standard Model (SM), such as supersymmetry and universal extra dimensions [1, 2, 3, 

j> ! 4, 5, 6]. An important uncertainty affecting the accuracy of any searches is the degree 
to which the detector interactions of hadronic SMPs can be modelled. In this paper a 



X 

^ model is presented for the scattering in matter of generic heavy hadrons containing either 
■ — ' up-like or down-like exotic quarks with charges ±|e and ±|e, respectively. The model is 
based on Regge phenomenology [7] and the Quark Gluon String Model (QGSM) [8]. 

Collider searches have already ruled out the direct pair production of a range of SMPs 
of masses up to around 200 GeV [1]. The LHC will open a new discovery window for 
SMPs with masses up to several TeV. The ATLAS and CMS experiments have already 
developed early SMP search strategies (see, for example, Refs. [9, 10, 11]), which usually 
require a track associated with a slow penetrating particle. To estimate the efficiency of 
such a search requires an understanding of the processes by which a SMP will interact 
with detector material. In the case of exotic heavy leptons, only electromagnetic energy 
loss could be expected. Hadronic SMPs could, however, also interact strongly, leading to 
additional energy loss and larger rates of SMPs stopped in detector material. Further- 
more, hadronic SMPs could undergo charge exchange reactions, which could lead to event 



^The term stable is taken to mean that the particle will not decay during its traversal of a detector. 
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topologies in which the exotic hadron appears to possess varying values of electric charge 
during its passage through the detector. Models of the interactions which a SMP may 
undergo are thus necessary in order to devise effective search strategies which can not 
only suppress SM background but also allow SMPs possessing different quantum numbers 
to be experimentally distinguished. 

While the electromagnetic interactions of massive objects are well understood [12], it 
is uncertain how hadronic interactions of an exotic hadron should be treated. Several 
models of varying sophistication have been proposed for the scattering of hadrons which 
comprise an exotic colour octet, eg gluino [13, 14, 15, 16]. There is so far only one detailed 
model [15, 16] of energy loss and charge exchange associated with the scattering of heavy 
hadrons containing exotic quarks. This approach, implemented within Geant [17], used 
a black disk approximation to obtain the total cross section and phase space arguments to 
predict the different types of reactions. By using Regge phenomenology and the QGSM 
we present in this work a complementary model which can also be included in Geant to 
aid future searches. 

This paper is organised as follows. The mass spectra of hadrons containing quarks 
and generic features of their scattering processes in matter are discussed. The QGSM and 
Regge-phenomenology are then used to provide a parameterisation in terms of relevant 
kinematic variables of single inclusive particle production arising from exotic hadron- 
nucleon scattering. This is used to estimate the average energy loss in such collisions and 
the expected rates of charge and baryon exchange processes. Using an ansatz of stable 
fourth generation quarks, we then describe the impact of this work on the prospects for 
detecting exotic hadrons at the LHC. 

2 Properties of exotic hadrons 

When considering interactions with matter, it is important to know the mass hierarchy 
of exotic heavy hadrons (hereafter referred to as i7-hadrons). This determines the states 
to which a i^-hadron, produced either in the primary interaction or after scattering with 
matter, would rapidly decay. Here, we consider only the lowest lying hadronic states, 
formed with an exotic heavy colour triplet charged object Q and light u and d quarks^. 

As outlined in [3, 15], the lowest lying neutral and charged mesonic states should be 
stable since the mass difference between them is expected to be far smaller than the pion 
mass. 

The baryon mass spectra, however, are more complicated [15]. The baryon state Qud 
containing light diquark system of spin is the lowest lying state HqucI with the other 
states Hquu and HQdd forbidden due to the requirement of anti-symmetric baryon wave 
function. The heavier states with diquark spin 1 {Hqud, Hq^d, and Hquu) are, however, 
all possible. In an analagous way to the charm sector in which the decays strongly 

to the Ac baryon [12], it would thus be expected that any produced heavy baryons would 
decay to the low mass state Hq^d via pion emission. This state would thus be charged 
(neutral) for i7-hadrons consisting of up-like (down-like) exotic quarks. As is shown in 

^Unless stated otherwise, charge conjugate states are also implied throughout this paper. 
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Section 4, this difference in charge of the stable if-baryon has a large impact on the 
observable rates of i?-hadrons at the LHC. 

In advance of a discovery it is premature to attempt to calculate precisely the rate and 

mass of each species of exotic hadron. Here, it is assumed that iJ-hadrons are degenerate 
in mass and that a sample of stable if-hadrons formed in high energy collisions, such as at 
the LHC, comprise 90% mesons, divided equally between charged and uncharged states, 
and 10% baryons. 

3 Interactions of i/-Hadrons in Matter 

Although phenomenological models of the type presented in this work are needed to 
predict some of the fine details of exotic hadron interactions in matter, it is nevertheless 
possible to build up a qualitative picture of the scattering process [15, 18]. Owing to the 
size of its wavelength, the heavy exotic quark will be a spectator, and it is the low energy 
hght quark system which interacts. Thus Regge phenomenology and the QGSM [8, 19] 
can be employed to model the interactions of exotic hadrons in matter. 

The QGSM is based on the 1/Nc expansion in QCD [20, 21] and partonic interpreta- 
tion of reggeon diagrams. All QCD diagrams are classified according to their topology. In 
this approach one can distinguish two classes of scattering processes: reactions mediated 
by (a) reggeon and corresponding to planar QCD diagrams and (b) pomeron exchange 
related to the cylinder-type diagrams in elastic scattering. Exotic hadrons containing a 
light constituent anti-quark, such as a Hq^ or a Hq^^ interact via pomeron and reggeon 
exchanges, the latter processes being due to the annihilation of hght antiquarks with 
the quarks of detector matter. Conversely, hadrons containing a light constituent quark 
(i^Qg, Hqqq ) can only interact via pomeron exchange. Multiparticle production related 
to reggeon and pomeron exchange processes are shown in Fig. 1. 

H ^ H' 




Figure 1: Diagrams representing reggeon (left) and pomeron (right) exchange between a 
if-hadron and a proton. 
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3.1 Total Cross Section of iif-hadron Scattering 

Let us consider the process of an interaction of a if-hadron with a nucleon of the target 
nucleus in the target rest frame. In this frame the hght antiquark of a if-hadron carries 
only a small fraction of the total energy E 

Eg ~ = ^mg± (1) 

where 7 = E/Mh and mq± is the transverse mass of the light antiquark. It was shown in 
the framework of the QGSM [8], that the planar diagram contribution to the total cross 
section cr^(s) is universal for the same energy of the annihilating antiquark. This implies 
that the contribution to the total reggeon cross section {aR{E)) can be written as: 

aniE) = Kapi{E = 7m,x) = KgR{2-fmg^/Eor^^'^-' (2) 

where K is the number of possible planar diagrams and Eq = 1 GcV. The vertex parameter 
gpt can be evaluated from the data on cross sections of hadronic interactions [22], and 
the intercept of the exchange degenerate Regge tragectories aji{0) is equal to 0.5. 

The pomeron contribution to the total cross section can be estimated with two models. 
In the additive quark model the P-contribution for the exotic meson is two times smaller 
than the corresponding value for pion-nucleon scattering since only one light quark (an- 
tiquark) is present. In the models in which hadrons are considered as colour dipoles the 
pomeron cross section is determined by the square of sizes of colliding hadrons. The mean 
radius squared of a hadron containing one heavy and one light quark is about 1/2 of 
of a hadron made of two light quarks and we get the same estimate for P-contribution 
as for the additive quark model. The pomeron cross section (ap) depends on energy as 

(TP ~ {2^m,jEor^^'^-' (3) 

The reggeon contribution to the cross section for a iJ-meson and a nucleon within a 
nucleus, which consists of equal amounts of protons and neutrons, can be derived as the 
difference between the reggeon contributions to a{7T~p) and a{7T~^p) data [7] multiplied 
with a factor 1.5. 

Fig. 2 shows the expected cross section for an exotic meson scattering off a stationary 
nucleon in a nucleus comprising equal amounts of protons and neutrons as a function of 
the Lorentz factor 7 of the exotic hadron. The contributions from reggeon and pomeron 
exchange processes are shown. 

Anti-baryons and baryons may interact via both reggeon and pomeron exchange, and 
pomeron exchange only processes, respectively. To obtain the overall cross sections for 
interactions involving baryons and anti-baryons, the pomeron contribution to the meson 
cross sections shown in Fig. 2 is doubled to take into account the extra light quark 
contribution. The reggeon contribution to anti-baryon scattering is taken to be twice the 
value for meson scattering with an added contribution from processes in which exotic anti- 
baryons can annihilate to exotic mesons and ordinary mesons. In the QGSM this process 
is described by planar diagrams with the annihilation of diquarks and string-junctions 
[23]. At large gamma the last process of string-junction annihilation dominates. In the 
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Figure 2: Pomcron (dashed) and reggeon (dotted) contributions to the exotic-meson- 
nucleon cross section. The sum of the two processes is shown as a sohd hne. 

same way as for usual annihilation it will decrease with energy as ~ 1/7^. At 7 ~ 1 the 
annihilation cross section can be large (~ 30 mb), and the annihilation cross section is 
thus taken to be 307~2 mb here. 

An exotic baryon can convert into an exotic meson in one of the sheets of the cylinder 
diagram of Fig. 2. This process can however be suppressed at energies close to threshold 
[10] by phase space effects and the absence of available pions within the nucleus which 
would be needed for reactions in which exotic baryons become mesons. 

The Regge approach is valid for 7 ^ 1, however it is known from experience with 
hadronic interactions that the Regge description works reasonably well on average for 
values of 7 ~ 1 [7]. The i7-hadron nucleon scattering cross section {(Thn) is used to 
estimate the cross section for the interaction of a if-hadron with a nucleus of atomic 
number A through anA = 1.25(jHAr^°'^ [17]. 

Fig. 3 shows the predicted distribution of the Lorentz factor for fourth generation 
quarks pair produced at the LHC, as calculated by Pythia [24]. Here is the total 
number of i^f-hadrons for each mass point: 200, 500, and 1000 GeV. As can be seen, for 
increasing mass the quarks will typically be produced with progressively smaller speeds. 
It should also be noted that a produced if-hadron must possess a speed above a certain 
threshold (for ATLAS this is conservatively estimated as 7>;1.4 [10]) to satisfy timing 
requirements for the particle to be triggered, read-out, and assigned to the correct bunch 
crossing. 

3.2 Differential Cross Sections and Energy losses 

In determining the kinematics of the scattering process, we consider the inclusive process 
H + N ^ H' + X , where H, H' , N and X are the incoming exotic hadron, the outgoing 
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Figure 3: The distribution of the Lorentz factor 7 of fourth generation quarks pair pro- 
duced at the LHC, as predicted by Pythia. The spectra are shown for quark masses: 
200, 500 and 1000 GeV. 



exotic hadron, the target nucleon, and whatever else is produced in the interaction, re- 
spectively. The kinematics of such an interaction can be specified by three independent 
kinematic variables. Commonly used variables are t, the usual four-momentum transfer 
between the incoming and outgoing exotic hadrons, s, the squared center-of-mass energy 
of the interaction, and Mx, the mass of the final state X. 

The final iJ'-hadron carries a fraction of energy xp close to unity and only a small 
fraction of energy 1 — xp ^ rnq^/Mn -C 1 is transferred to production of hadrons. This 
justifies the application of the triple-regge formulae to provide a description of inclusive 
cross sections. Strictly speaking the triple-regge description is valid for m\ ^ IGeV"^ 
and the rapidity difference between H' and rest hadrons Ay > 1. This is equivalent to 
the condition 2'ymq±mi\f /Mj^ ^ 1. In hadronic interactions, the triple-regge description 
works usually up to Ay ~ 1 and we will assume in the following that the same is true 
for interaction of if-hadrons [7]. A representation of rapidity gaps between the quark 
systems in the three reggeon exchange diagram is shown in Fig. 4. 

Expressions for the contributions of different triple-regge terms iik (Fig. 4) to inclusive 
cross sections is straightforward to obtain, noting that for reggeons i corresponds the factor 
exp(2(«j(t) — l)Ay), while an exchange by the reggeon k leads to the factor exp((Q;fe(0) — 
l)yq)- Here, yg = ln{Mjr/{mq±mN) is the rapidity interval covered by produced hadrons 
(the total rapidity Y — \n(2E /Mh) — ln(27) = Ay -\- yq). As for the total cross section 
we consider as exchanged reggeons i, k the pomeron P and secondary reggeons R. Thus 
we have the following triple-regge contributions: RRR,RRP,PPR and PPP. 

Using the rules described above we can write inclusive cross sections for the corre- 
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Figure 4: Left: triple-regge diagram describing the process H + N ^ H' + X. Right: 
representation of rapidity gaps between the quark systems in the three reggeon exchange 
diagram. 



sponding triple-regge terms in the following forms: 



dtdMx 



dtdMj. 



2 (7, Mi) 
(7, Ml) 



dtdMj. 



X, 



1 21 M"^ f M'^\ 

^4(7)CiJ«iiexp[(2B«^ + Burr + 2a'^\n{^^))t] ^ 

j^ct'MCrrp exp[{2BRH + Brrr + 2a'p 
j^aU^)CppRexp[{2BpH + Brrr + 2a'p 

X X \ X / 

j^aU^)CpppeM{'^BpH + Bppp + 2a'pH^))t] (7) 

X X \ X / 



(4) 
(5) 
(6) 



where = ^^(0) - 1= -0.5, Ap = ap(0) - 1= 0.12, a'^ 



R 



0.9 Ge1/-2, a'p 



0.25 

[8] and = m^vm^x = 0.5 GeV^. 

The parameters, Cuj, and Baj can be determined using Regge factorization from the 
triple-regge description of inclusive spectra in high-energy hadronic interactions [19]. Let 
us emphasise that the RRR-term corresponds to the diagram of Fig. 1 (left), which repre- 
sents the cutting of the planar diagram or R-exchange, while the RRP-term corresponds 
to the cutting of the cylinder-type diagram in Fig. 1 (right). Due to conservation of H- 
hadrons integrals over and t give ap and ap contributions to the total cross section 
correspondingly. 

The PPR and PPP-terms describe the diffractive dissociation of a nucleon and their 
cross sections can be calculated, using factorization from the corresponding cross sections 
extracted from pp-inter act ions 



a 



ppi _ (^p{Hpf pp, 



Hp 



apippf 



a. 



pp 



7 



Here, we neglected the small difference in t-dependence for Hp and PP vertices. Taking 
into account that "'^^f^^^ pa 1/4 and that the sum of PPR and PPP-contributions for 
pp-coUisions in the relevant energy domain does not exceed 2mb, we obtain very small 
cross sections for diffraction dissociation of a nucleon in ifp-interactions: 0.12 mb. Thus 
these cross sections constitute only about 1% of the total cross section and can be safely 
neglected in our estimates of energy losses. 

For parameters characterising the t-dependence of RRR and RRP-terms we take the 
same values as have been extracted from analysis of pp-interactions [19]: 2BRH + BjiRji — 
2Brh + Brrp = 4 GeV-^ 

The energy loss of a i7-hadron is given by: 

AS = ^^-"■" + 1" (9) 
The average energy loss can thus be calculated: 

/ rp\ _ ■Imff+mTr ^ •'\t\min I I d\t\dMx ('XtW 
rMxmax J/l/f r\t\max ^Ul iPa 
JmN+mn^-^^^X J\t\min "l''ld|t|dMx 

Here, mjv and ttIt^ were taken as the mass of the proton and a charge pion, respectively. 

The upper limit on Mx is taken to be the lower of the following two limits: Mxmax = 
(27Mq)2, which represents the condition Ay = or Mxmax = y/s — mn from energy- 
momentum conservation and where ttih is the mass of the interacting if-hadron. 

The limits on t are given by 

\t\min,max{Mx) = 2[E{mN)E{Mx) T p{m^)p{Mx) - mW (11) 
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where E{m) = p{m) = ^^(^'^Ji'^'\ and A(a, 6, c) = + + - 2{ah + 

ac + he) . 

Fig. 5 shows the mean energy loss associated with two triple-regge contributions (RRR, 
RRP) as a function of 7. The two contributions give similar values except at the very 
highest values of 7, at which the production cross section for iJ-hadrons is tiny. It was 
also found that the energy loss at a fixed jS is insensitive to the i7-hadron mass for mass 
values greater than around 10 GeV. The principal uncertainty in the total energy loss 
comes from the uncertainty in the value of m^^, which enters via the parameter Mq 
and was chosen to be 0.5 GeV^ for this work. Taking values of Mq of 0.3 or 0.7 GeV^ 
changes the average energy values typically by around 50% and this would not change 
the conclusions of this paper. 

For calculations of energy losses of a given i7-hadron type we need in principle to 
know only the losses associated with each triple-regge term and the relative weights of 
different terms. For RRR and RRP-terms these weights arc given by cr^(7)/o"*''*(7) and 
o"p(7)/(T*°*(7) correspondingly. However, since the difference in energy loss from the two 
contributions is anyway very small, iif-hadron energy loss for both reggeon and pomeron 
interactions was assumed to be that given by the RRP term. 

A Monte Carlo method was used to extract the predicted distributions f{E) — 
of energy loss {E) at a range of values of (3, where N is the total number of i?-hadrons 
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Figure 5: Average energy loss per interaction as a function of 7 for two triple-regge terms. 



at each /3 point. The functional form f{E) for RRR and RRP terms in the if-hadron 7 
range expected at the LHC (7;$ 10) can be approximated by Equation 12. 

f{E) = exp[-10(S - 0.3)^(1 - pn ^^^^o_^y, (12) 

This function was used to arrive at the results presented in Section 4. Fig. 6 shows f{E) 
for an exotic i7-hadron of mass 200 GeV for four values of /3. The distribution typically 
extends up to around a few GeV, though it is peaked at around 0.5 GeV. 

3.3 Charge Exchange and Baryon Formation 

Charge exchange processes will naturally take place via the formation of light quark pairs 
from the vacuum both for planar (Fig.l) and cylinder-type (Fig. 2) diagrams. As u and 
d quarks arc produced with equal probability we expect that charge exchange happens 
with 50% probability per interaction. 

Within the Lund string model [25] exotic baryons only account for around 10% of all 
exotic hadrons produced at the primary interaction, a prediction which is also made in the 
QGSM [26] . Baryons can also be produced as a result of hadronic interactions of exotic 
mesons in matter. The probability that a given inelastic collision involving an exotic 
meson results in baryon formation takes place is taken to be 10%, a value motivated by 
investigation of baryon production in QGSM and low energy hadron scattering data [8]. 
The proportion of different baryon species formed in such interactions is taken to be the 
same as at the primary interaction. 
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Figure 6: The predicted distribution of the energy loss per interaction of a i7-hadron 
of mass 200 GeV and for different values of p. 



4 Experimental Signatures 

In this section the model of hadronic scattering outlined in Section 3, together with the 
well established theory of electromagnetic energy loss of charged particles in matter [12], 
is used to estimate the energy losses and rate of charge exchange reactions of a iJ-hadron 
propagating through the ATLAS detector [27] at the LHC. Calculations of the cross section 
for fourth generation quark pair production at the LHC are then used to predict the rates 
of various event topologies associated with i7-hadron production. The event topologies are 
characterised by the values of charges possessed by the if-hadrons in the inner and muon 
tracking chambers, i.e. before and after scattering in the ATLAS calorimeter system. 

4.1 Energy Loss and Charge Exchange 

To estimate the kinematic distributions of if-hadrons produced at the LHC, the Pythia [24] 
program was used to generate samples of 50,000 fourth generation quark pair production 

events for quark masses 200, 500 and 1000 GeV. It was assumed that the proportion 
of different types of stable if-hadrons formed from these quarks follows the prescription 
given in Section 2. To ensure that if-hadrons would belong to a high acceptance re- 
gion of the ATLAS detector, the initial value of P of the iJ-hadrons was required to be 
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greater than 0.7 and the pseudorapidity was restricted to |?7| < 2.5 [10]. The remain- 
ing if-hadrons were then transported, using a Monte Carlo method, through material of 
thickness corresponding to the part of the ATLAS detector enclosed by the muon tracking 

chamber. 

The material distribution of the sub-systems of the ATLAS detector enclosed by the 
muon detector varies between around 11 and 19 interaction lengths as a function of pseu- 
dorapidity [27]. The largest material contribution (~ 70% of the total thickness) arises 
from the absorbing material in the different hadronic calorimeter systems. These systems 
arc the Tile Calorimeter (TileCal) [28] and Hadronic EndCap Calorimeter (HEC) [29], 
which use iron and copper, respectively, as absorbers. Thus, for a generated i/-hadron 
falling in the pseudorapidity regions |?7| < 1.5 (1.5 < |?7| < 2.5) covered by the TileCal 
(HEC), iron (copper) was used to represent the ATLAS detector material in the Monte 
Carlo calculations presented here. 

Fig. 7 shows differential distributions related to the energy loss and interactions of 
if-hadrons as they pass through the detector material. Spectra are shown separately 
for i^-hadrons formed from different types of exotic quarks and anti-quarks with masses 
200 and 1000 GeV. The distributions are normahsed to the total number N oi a, given 
type of if-hadron satisfying the /3 and rj requirements. In the following discussion of 
the plots Hq[Hq) is used as a generic term to denote a if-hadron with an exotic quark 
(anti-quark) while the terms Hu, H^, Hfj, and Hjj are used to denote i?-hadrons with 
up-like and down-like exotic quarks and anti-quarks. 

The distributions of total energy loss are peaked at ~4-5 GeV and extend up to around 
10 GeV for the different iJ-hadron types, with no substantial mass dependence. The Hd- 
hadrons typically show lower energies than the i^iy-hadrons. This is due to differences in 
ionisation energy loss since ifjj-hadrons are more likely to propagate through the material 
with zero electric charge. A if^-hadron can, for example, start as a neutral meson and 
then be converted into a neutral baryon. The effect can be seen as a peak at low ionisation 
energy low for if£)-hadrons. Smaller peaks at low ionisation energy loss are also visible 
for the i^Q-hadrons, which arise from events in which those i^-hadrons propagate through 
the material mostly as neutral mesons. 

The hadronic energy loss for the if-hadrons decreases slightly with mass owing to the 
typically lower speeds of the more massive iJ-hadrons (see Fig. 5). The hadronic energy 
loss for i^Q-hadrons peaks at around 3-4 GeV, which is larger than that for iJ^-hadrons 
(~ 1 GeV). This can be understood as a consequence of the different reaction channels, 
which are open for Hq and iJg-hadrons. iJ-hadrons are dominantly mesonic following 
hadronisation and, although they can form baryons, still interact mostly as mesons. As 
explained in Section 3, Hg-mesons may scatter via reggeon or pomeron exchange, unlike 
i^Q-mesons, which interact with a lower cross-section since only pomeron exchange is 
possible in this case. This is seen in the distribution of the multiplicity of hadronic 
interactions, which peaks at ~ 5 (~ 2) for iJ-hadrons with exotic quarks (anti-quarks). 

Within the if-hadron kinematic region under study, the stopping of if-hadrons is 
negligible, as can be understood from their low energy losses. The decrease in speed Sp 
of a i?-hadron following the traversal of detector material in the ATLAS detector was 
calculated and it was found that 5/3<0.005 (5/3<0.02) for a mass of 200 GeV (1000 GeV). 

The results given above were obtained assuming that a neutral i^-hadrons could not 
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Figure 7: Distributions of energy loss and hadronic scattering for iJ-hadrons of masses 
200 and 1000 GeV and for exotic quarks of charges ±|e and ±|e. The left (right) column 
represents /7-hadrons containing an exotic quark (anti-quark). Distributions of the total, 
lonisation and hadronic energy loss is shown along with the multiplicity of interactions. 
The distributions assume no mixing of neutral i/-mesons. 
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oscillate into its anti-particle, something which has been studied in the context of su- 
persymmetry [3, 30]. In a strictly generic scenario of new physics, as considered here, 
it is prudent to consider the possibility of two extreme cases: maximal and no mixing. 
Allowing maximal mixing produces distributions (not shown), which are very similar to 
those given in Fig. 7. However, together with the effects of hadronisation and scattering 
interactions it can give rise to striking event topologies, which can be used to detect and 
characterise if-hadron, as is shown in Section 4.2. 

4.2 Expected Rates of if-hadrons 

Fig. 8 shows the expected cross section for the pair production of fourth generation quarks 
of a single flavour at the LHC, as predicted by Pythia. The cross section is seen to 
fall steeply with mass. Nevertheless, for a full year's LHC running at low luminosity, 
corresponding to 10 fb~^, around 25,000 pairs would still be expected for a 500 GeV 
mass, implying that exotic hadrons can be sought with early LHC data. 
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Figure 8: Predicted cross section, calculated to leading order, for the pair production of 
fourth generation quarks at the LHC. 

Following the acceptance cuts and iJ-hadron transport through ATLAS described in 
Section 4.1, Tabs. 1 and 2 contain the expected yields for a range of event topologies 
arising from the pair production of exotic quarks possessing ±|e and ±|e. Rates are 
shown for if-hadrons of mass 200, 500 and 1000 GeV and an integrated luminosity of 
10 fb~^ is assumed. The topologies are: 

1. At least one i/-hadron is produced with non-zero charge and retains that charge 
value throughout its passage in the detector material implying that the inner and 
muon tracking chambers measure the same charge; 
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2. Both particles are produced with non-zero charge and retain those charge values. 

3. Both particles are produced with non-zero charge and one retains its non-zero charge 
but the other is converted into a neutral state. 

4. One particle is produced with zero charge whilst the other has non-zero charge; the 
particle produced with zero charge converts into a charged state whilst the other 
retains its original value of charge. 

5. Both particles are produced with non-zero charge but leave the detector material as 
neutral states. 

6. At least one particle leaves the detector material possessing a non-zero electric charge 
of opposite sign to the charge with which it was produced. 





No mixing 




Maximal mixing 


Topology 


Mass (GcV) 




Mass (GeV) 






200 


500 


1000 


200 


500 


1000 


1 


4.9 X 10^ 


4.3 X 10=^ 


57 


4.1 X 10^ 


3.5 X 10=^ 


48 


2 


3.0 X 10^ 


2.6 X 10^ 


3 


2.2 X 10^ 


1.9 X 10^ 


2 


3 


9.6 X 10^ 


8.3 X 10^ 


9 


8.2 X 10^ 


6.8 X 10^ 


8 


4 


6.0 X 10^ 


5.2 X 10^ 


6 


4.8 X 10^ 


4.0 X 10^ 


5 


5 


6.4 X 10^ 


5.3 X 10^ 


6 


6.3 X 10^ 


5.5 X 10^ 


6 


6 











8.1 X 10^ 


7.2 X 10^ 


9 



Table 1: Expected rates of various topologies, corresponding to an integrated luminosity 
of 10fb~^, for different event topologies arising from the pair-production of exotic quarks 
of charge ±|e. 

As can be seen, substantial rates for the mass points 200 and 500 GcV can be expected 
for each of the topologies. While the first two topologies offer the classic signature of a 
single or a pair of high momentum penetrating particles, scenarios 3 and 4 show evidence 
of charge exchange and can thus be used in combination with calorimeter information to 
identify the SMP as being hadronic and reject alternative scenarios, such as stable Icptons. 
Scenario 5 may be the most experimentally challenging topology if a muon-bascd trigger 
is used, although, as seen, this topology accounts for only a small fraction of events. The 
observation of topology 6 could demonstrate mixing in the if-hadron sector. 

A comparison of Tabs. 1 and 2 shows the rates for all of the considered topologies bar 
scenario 5 are greater for iJ-hadrons with up-like exotic quarks than those for those with 
down-like exotic quarks. However, the situation is reversed for topologies (5). This can 
be explained as being due to the production of baryons which have zero (non-zero) charge 
for ±|e (±|e) charge exotic quarks. 
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No mixing 




Maximal mixing 




TonoloETv 


Mass (GeV) 




Mass (GeV) 






200 


500 


1000 


200 


500 


1000 


1 


8.8 X 10^ 


8.0 X 10^ 


107 


8.0 X 10^ 


7.2 X 10^ 


97 


2 


1.2 X 10^ 


1.0 X 10^ 


12 


8.6 X 10^ 


7.7 X 10^ 


8 


3 


1.3 X 10^ 


1.2 X 10^ 


13 


1.0 X 10^ 


8.9 X 10^ 


10 


4 


1.7 X 10^ 


1.5 X 10^ 


17 


1.8 X 10^ 


1.5 X lO-'^ 


17 


5 


2.9 X 10^ 


2.5 X 102 


3 


2.5 X lO'^ 


2.1 X 102 


2 


6 











1.9 X 10^ 


1.7 X 10^ 


23 



Tabic 2: Expected rates of various topologies, corresponding to an integrated luminosity 
of lOfb^^, for different event topologies arising from the pair-production of exotic quarks 
of charge ±|e. 

4.3 Detector Effects 

It is beyond the scope of this paper to study in detail the expected response of the ATLAS 
detector. Such work is more appropriately made with the ATLAS simulation and event 
reconstruction packages [31]. However, a few comments can be made based on earlier work 
employing such programs and the results presented in this Sections 4.1 and 4.2 should be 
considered in the light of these. 

It is desirable for any search to include calorimeter information. The observation of 
small hadronic energy depositions associated with a penetrating SMP could be used, for 
example, to reject backgrounds of muons. Furthermore, as shown in Fig 7, differences 
in hadronic energy loss between a pair of SMP candidates could provide evidence for 
stable exotic quarks. To study the feasibility of exploiting calorimctry in this way, the 
model presented here should be implemented within a package such as Geant in order 
to include effects such as nuclear fragmentation [17]. 

A conservative estimate of the trigger efficiency for a slow moving particle using a 
muon trigger of 50% has been made [10]. Therefore, with the exception of topology 5, 
events corresponding to the topologies listed above should be recorded in this way. It 
may be possible to select events corresponding to topology 5, through a jet or missing 
transverse energy trigger, in the case when one of exotic quarks emits a high transverse 
momentum quark or gluon. The track reconstruction efficiency in the inner and muon 
tracking chambers should be over 90% in the kinematic region under investigation [27], 
though further losses could occur due to charge exchange interactions of iJ-hadrons in the 
different tracking systems. More studies on efficiency losses using, for example, the model 
presented in this work is desirable. Furthermore, the charge misidentification probability 
for tracks in the muon (inner) detector should be around 2% (4%) for transverse momen- 
tum pt values up to around 1 TeV and falls quickly with decreasing pt [27]. While it could 
be expected that the charge misidentification probability may be degraded for interacting 
i7-hadrons, earlier work with gluinos involving a full ATLAS detector simulation showed 
that a signature of an exotic hadron apparently reversing the sign of its charge is a useful 
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search observable [32]. We would thus expect topology 6 to be distinguished. 

SM backgrounds can be suppressed in a number of ways, for example by making 
requirements on the SMP's transverse momentum {pt^IOO GeV) and event shape vari- 
ables [10, 32]. The most promising method is to use a time-of-fiight technique [10] with 
which exotic hadron searches in early LHC data may only suffer from a handful of back- 
ground events. 

5 Conclusions 

Heavy exotic quarks are predicted in a number of scenarios of physics beyond the SM. 
This paper presented a generic model based on Regge theory and the Quark Gluon String 
Model to describe the interactions of heavy hadrons in matter. The work showed how 
strong interactions of exotic hadrons may provide useful observable to detect and quantify 
the properties of any exotic stable massive particle which may be observed. Distributions 
relating to the interactions and estimates of expected event topologies at the LHC were 
shown. A natural next step in this work would be the implementation of this model within 
Geant and its subsequent use in searches at colliders and studies of discovery potential. 
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